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Abstract 

A search for ^ A^A*(892)^ decays is performed using pp collision data, 

corresponding to an integrated luminosity of 1.0 fb~^, collected with the LHCb 
detector at a centre-of-mass energy of 7 TeV. The —)• ^^^*(892)*^ decay is 

observed for the first time, with a significance of 7.1 standard deviations. The 
branching fraction is measured to be 

^ K°K*{892f) + ^ K°K*{892f) = (16.4 ± 3.4 ± 2.3) x 10“®, 

where the first uncertainty is statistical and the second is systematic. No evidence is 
found for the decay B^ —?• A^A*(892)^ and an upper limit is set on the branching 
fraction, B{B° K^K*{892f) + B{B^ K^K*{892f) < 0.96 x 10"®, at 90% 
confidence level. All results are consistent with Standard Model predictions. 
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1 Introduction 


Violation of the combined charge-conjugation and parity symmetry (CP) is one of the fun¬ 
damental ingredients to explain a dynamical generation of the observed matter-antimatter 
asymmetry in the universe |^. In the Standard Model of particle physics (SM), CP 
violation in the quark sector is generated by a single complex phase in the Cabibbo- 
Kobayashi-Maskawa matrix [^[^. However, the observed baryon asymmetry is too large 
to be explained by the SM mechanism alone [^. Non-leptonic B meson decays dominated 
by amplitudes involving a quark and a W boson in a loop (penguin) are sensitive to the 
presence of non-SM physics processes. These processes could provide additional sources of 
CP violation that might explain the observed baryon asymmetry. The —?• K^h^h ^ 

(h, h' = TT, K) decays are interesting for CP violation measurement^ Knowledge of the 
branching fractions of the various sub-modes, as reported in this paper, is an important 
input to the theory of CP-violation, particularly models of new-physics contributions 
to 6 —)• s transitions Ref. [^| . The measurements also allow tests of QCD models (see, 
for example, the predictions in Refs. [6]-[^). If sufficient data are available, a common 
approach for three-body decays is to perform an amplitude analysis by studying the 
structure of the Dalitz plot |^. If data are less abundant and the decay products originate 
from intermediate resonances, as in the present analysis, a quasi two-body approach can 
be used. 

The LHCb collaboration has provided results for inclusive —)■ K^h^h ^ decays 

and more recently the hrst measurements of meson decays to K*(892)“7r 
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and 


K*(892) hnal states d- An initial search for the neutral decay B^ —)■ K* (892 )^has 


12 


. In this paper a search for Kg K* (892)'- 

meson, hereafter denoted by decays to the 
invariant mass region around 


been reported by the BaBar experiment 
decays is reported, where the K*(892)° 

K^ir~ final state. The resonant structure in the K+tt” 

1 GeV/C is analysed to determine the number of decays that proceed through an interme¬ 
diate K*^ resonance. The branching fraction is measured relative to the B^ —)■ Kg 7r’''7r“ 
decay [T^ , using the relation 


A'»ir+ir-) 


€ 30 - 


(g) S ^ _ 

S (s. 
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^K^K*o 
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fdis) ' H(K*o^K+7r-)’ 
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where N represents the number of observed decays, e the total efficiency, and fs/fd the 
ratio of the fragmentation fractions of a 6 quark into a or a meson 14 ^ 

B{K*^ —)■ K’''7r“) is the branching fraction of the K*° meson into K+vr” final state, 
the following, the —)■ KgK*° and B^ —)■ decays are referred to as signal and 

normalisation channels, respectively. 


and 

In 


^Charge-conjugate modes are implicitly included throughout this paper. 
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2 Detector and simulation 


The analysis is performed using pp collision data recorded with the LHCb detector, 
corresponding to an integrated luminosity of f.Ofb”^, at a centre-of-mass energy of 
7TeV. The LHCb detector f7, f8 is a single-arm forward spectrometer covering the 


pseudorapidity range 2 < p < h, designed for the study of particles containing 6 or c 
quarks. The detector includes a high-precision tracking system consisting of a silicon-strip 
vertex detector surrounding the pp interaction region, a large-area silicon-strip detector 
located upstream of a dipole magnet with a bending power of about 4Tm, and three 
stations of silicon-strip detectors and straw drift tubes placed downstream of the magnet. 
The tracking system provides a measurement of momentum, p, of charged particles with a 
relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200GeV/c. The 
minimum distance of a track to a primary vertex (PV), the impact parameter, is measured 
with a resolution of (15 -|- 29/pt) M-ni, where px is the component of the momentum 
transverse to the beam, in GeV/c. Different types of charged hadrons are distinguished 
using information from two ring-imaging Cherenkov detectors. Photons, electrons and 
hadrons are identihed by a calorimeter system consisting of scintillating-pad and preshower 
detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identihed 
by a system composed of alternating layers of iron and multiwire proportional chambers. 

Simulated events are used to determine the efficiency of the selection requirements, to 
study possible sources of background and to determine the parametrisations used to model 
the data. In the simulation, pp collisions are generated using Pythia 6 19 with a specihc 
LHCb conhguration 20 . Decays of hadronic particles are described by EvtGen 


in which hnal-state radiation is generated using Photos [^. The interaction of the 
generated particles with the detector, and its response, are implemented using the Geant4 
toolkit 


23 as described in Ref. 24 


3 Event selection 

The online event selection system (trigger) consists of a hardware stage, based on 
information from the calorimeter and muon systems, followed by a software stage, in which 
all charged particles with px > 500MeV/c are reconstructed. The hardware hadron trigger 
requires a calorimeter cluster with transverse energy greater than 3.5 GeV. In the offline 
selection, candidates are divided into two non mutually exclusive categories based on the 
hardware trigger decision. One category consists of candidates whose decay products 
satisfy the hadron trigger requirements, while the other consists of candidates from events 
in which other particles meet the hardware trigger requirements. Only events that fall into 
either of these categories are used in the subsequent analysis. The software trigger requires 
a two-, three- or four-particle secondary vertex with a significant displacement from the 
primary pp interaction points. At least one charged particle must have px > 1.7 GeV/c 
and be inconsistent with originating from any PV. A multivariate algorithm is used 
for the identihcation of secondary vertices consistent with the decay of a 5 hadron. 
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In the offline selection the —)■ K^K*^ decays are reconstrncted throngh the K*^ 

—^K^Tr~ and —j-tt+tt” decay modes, where the Kg candidate is constrained to its 

known mass and the B candidate is constrained to originate from a PV. Decays of 
Kg mesons are reconstrncted in two mntnally exclnsive categories: long Kg candidates, 
which decay snfhciently early that their danghter pions are reconstrncted in the vertex 
detector; and downstream Kg candidates, which have danghter particles that are only 
reconstrncted in the rest of the tracking system. As these two categories have different 
backgrounds, and the long Kg mesons have better momentum and vertex resolutions, the 
selection criteria for long and downstream Kg candidates differ. The selection criteria 
follow those in Ref. 
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Fully reconstructed background decays that have the same hnal state as the signal 
include contributions from B decays to hnal states involving charm mesons, such as Dh, 
with a Kgh~^h~ hnal state, or A° decays to A+h“, with A+ —)■ Kgp, where the proton is 
misidentihed as a or K^. In addition, B decays with an intermediate charmonium state 
like B^ —)■ J/ipKg, with , K~^K~, , can be present in the mass region of 

the normalisation channel. To reduce the contamination from these backgrounds, a veto is 
applied on the invariant mass of each of the possible intermediate states reconstructed 
under the corresponding hypothesis. Candidates are excluded if the reconstructed mass 
of a two-body intermediate state is within 30MeV/c^ (48MeV/c^) of the known mass of 
the relevant intermediate charm (charmonium) resonance of one of the backgrounds 
considered. No particle identihcation information is used at this stage. 

If a hnal-state hadron is misidentihed, signal yields can potentially be affected by decays 
into any Kgh^h'^ hnal state, especially when the h^h'^ proceeds through a resonance. 
Particle identihcation requirements on the two tracks originating from the B decay vertex 
are used to separate pions, kaons and protons, and to reduce this background to a negligible 
level. The largest source of background is due to random tracks that form candidate 
B or Kg decay vertices. A multivariate discriminant based on a boosted decision tree 
(BDT) algorithm 27,28 is used to reduce this background. The greatest discrimination 
in the BDT is provided by kinematic properties of the B meson, its hight direction 
with respect to the PV, and variables dehned analogously for its decay products. The 
optimisation of the BDT is described in Ref. 10 ; the selection requirement on the BDT 
response for this analysis is chosen to maximise e/(a/2 + ^/Nb) (^. Here, e is the signal 
efficiency, B represents the number of background events in the signal mass interval, which 
is estimated using data by extrapolating the number of background events from the upper 
mass sideband into the signal region, and a = 5 is the chosen target signal significance. 

The efficiencies are determined from simulation, except for the particle identihcation 
efficiencies. The latter are determined from data using samples of kinematically identihed 
charged particles from )■ D^7i~^ with )■ K~tt^, and A—)■ n~p decays, reweighted 
to match the kinematic properties of the signal. The BDT selection efficiency for signal 
is approximately 85% (90%) for downstream (long) signal decays; approximately 88% 
(95%) of backgrounds in the respective categories are rejected. The B^ —)■ Ag 7r’''7r“ decay 
selection efficiency is taken from Ref. 


10 . The efficiencies for the normalisation cha nn el 


are determined in bins of the Dalitz plane and are reweighted from data using the sPlot 
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4 Fit model 


Two-dimensional extended maximum likelihood fits to the unbinned and K^tt~ 

mass distributions are used to determine the event yields for the signal channel, while 
an independent one-dimensional ht to the mass distribution is used for the 

normalisation channel. The correlation between the two signal mass distributions is 
checked on simulation. The results do not show signihcant correlations and therefore 
the correlation terms are neglected in the ht. Candidates in the long and downstream 
categories are htted simultaneously. The signal ht is restricted to candidates in the mass 
regions 5000 < m{K^K^TT~) < 5800 MeV/c^ and 650 < m{K^TT~) < 1200 MeV/c^. The 
ht model consists of signal, non-resonant background, partially reconstructed background 
and combinatorial background components. 

The and the components of the signal are both parametrised as two Gaussian 
distributions with a power-law tail on each side. For each component the two functions 
share the peak position and the width parameters. The parameters describing the tails 
are determined by hts to simulated samples and subsequently hxed in the ht to data. The 
systematic uncertainty associated with this choice is found to be negligible. The B^-B^ 

The K*^ mass distribution is parametrised 


mass difference is hxed to the known value 13 


by a relativistic Breit-Wigner function with the peak position and width allowed to vary 
in the ht. 

The components in the B mass model that are non-resonant in K^tt~ are parametrised 
by the same function as the signal, sharing their peak positions and widths with the 
signal functions. The tail parameters are hxed according to the values obtained from 
simulation. The non-resonant component of the K^ti~ mass distribution is approximated 
by a normalised linear function as in Ref. 11 , with the zero point of the function on the 


abscissa determined by the ht. While the ratio between the non-resonant and the signal 
components is hxed to be the same for the two meson categories, it is independent for 
the B^ and the 5° candidates. 

Backgrounds from partially reconstructed decays are classihed into two categories. 
Decays such as R —)■ Dh are parametrised by means of ARGUS functions convolved with 
a Gaussian function in the B candidate mass, and linear functions in the K*^ candidate 
mass. The choice is based on simulation studies and previous hndings 10 ,^. In decays 
such as —)■ where one resonance decays as )■ iF+vr" while the other decays 

as K —)■ Kgir^, the B^ mass distribution is described using the same parametrisation as 
for the previous background, while the invariant mass distribution for K*^ candidates is 
described by a relativistic Breit-Wigner function sharing the peak position and widths 
with the signal component. The yield for these components are determined in the ht to 
data. 

The combinatorial background is modelled by an exponential function in the B can¬ 
didate mass distribution and a linear function in the K*^ candidate mass distribution. 
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Table 1: Signal yields obtained from the fits to and mass distributions and 

corresponding efficiencies. Only statistical contributions to the uncertainty are reported. 


Decay 

Downstream 


Long 

Yield 

Efficiency (%) 

Yield 

Efficiency (%) 

5°^ K^K*^ 

21 ± 6 

0.0174 ±0.0012 

25 ± 6 

0.0121 ±0.0008 

B°^ KlK*^ 

2± 3 

0.0183 ±0.0013 

1± 2 

0.0125 ±0.0009 

B°^ KItt+tt- 

828 ± 41 

0.0336 ±0.0010 

341 ± 23 

0.0117 ±0.0009 


These functions are found to give good agreement with the distributions in the appropriate 
data sidebands. The slopes of the exponential functions are independent for the long and 
downstream categories, while the abscissae of the linear functions are the same. All these 
parameters are allowed to vary in the fit. 

The parametrisation used to model the —)• AT®tt+tt” normalisation and the back¬ 

ground follow those used to fit the signal mode. In addition two other categories of partially 
reconstructed backgrounds are included; decays such as —)■ A'gTr+vr^y or B^ —)■ 

with r^—)• p° 7 ; and misidentified B^^^ —)• decays. Their parameters are fixed in 

the fit to the values derived from simulated samples. 

The observed K^K*^ and Kg7r^7r~ mass distributions and the corresponding fits are 
shown in Figs. and respectively. The signal yields are reported in Table [T] The 
B^ mode is dominated by the non-resonant component. The statistical significance of 
the Bg signal is determined using Wilks’ theorem 31 and by combining the long and 
the downstream samples. The significance including relevant systematic uncertainties, 
estimated by repeating the procedure with the signal likelihood convolved with a Gaussian 
function of width equal to the sum in quadrature of the systematic uncertainties, is 7.1 
standard deviations. 


5 Systematic uncertainties 

The model used to fit data and the limited knowledge in the efficiency determination 
are possible sources of systematic uncertainty. Many parameters in the fit are fixed to 
values obtained from fits to simulated data. The associated systematic uncertainties are 
determined from fits to pseudoexperiments generated assuming alternative values of the 
relevant parameters, corresponding to variations within uncertainties around their default 
values. The average difference between the yields determined in the pseudoexperiments 
and the nominal value is taken as a systematic uncertainty. 

The fit model does not account for the possible interference between the P wave of 
the K*{892)^ resonance and the S wave from other intermediate states, e.g. the non¬ 
resonant component or the 77* (1430)*^ resonance. The associated systematic uncertainty is 
determined by exploiting the distribution of 0k*o, defined as the angle between the flight 
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Figure 1: Distribution of (left) mass and (right) mass for signal candidates 

with ht results overlaid for (top) downstream and (bottom) long categories. The data are shown 
as black points with error bars. The overall fit is represented by the solid black line. The 
and signal components are the black short-dashed and dotted lines respectively, while 
the non-resonant components are the magenta short-dashed and dotted lines. The partially 
reconstructed backgrounds are the red triple-dotted line {B — Dh) and the blue triple-dotted 
line {Bg^ K*^K*^). The combinatorial background is the green long-dash dotted line. 


direction of the in the K*^ rest frame with respect to the direction of the boost from 
the laboratory frame to the K*^ rest frame. The cos distribution is described by a 
parabola, where the second-order term represents the signal P wave, the constant term is 
related to the S wave and the hrst-order term accounts for the interference. Using the sPlot 
technique , the cos 6k*o distribution of the signal P and S wave is unfolded from the 
other background components. A fit in the region of positive cos 9^*0 is performed using 
a second-order polynomial and the systematic uncertainty is determined as the relative 
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Figure 2: Distribution of mass for signal candidates with fit resnlts overlaid for (left) 

downstream and (right) long categories. The data are described by the black points with error 
bars. The overall fit is represented by the solid black line. The and signal components 
are the black short-dashed and dotted lines. The misidentified B^ decay is the black dashed line, 
respectively. The partially reconstructed backgrounds are the red triple-dotted line {B —)• Dh), 
the blue triple-dotted line {B^ —)• the violet dash single-dotted line {B^ —)• r]'K^) 

and the pink short-dash single dotted line {B^ —)■ it'g7r"’'7r“7). The combinatorial background is 
the green long-dash dotted line. Some of the contributions are small in the figures. 


difference between the integral of the function when the constant coefficient is allowed to 
vary and when this coefficient is fixed to zero. Due to the limited size of the )■ K^K*^ 
sample, the relative uncertainty obtained for the —>■ K^K*^ decays is also applied to 
the 5°—)- KgK*^ decays. 

Potential biases that may be associated with the maximum likelihood estimator are 
investigated using pseudoexperiments. The systematic uncertainty is determined as the 
average difference between the nominal value and the htted yields in the pseudoexperiments. 

The impact of the limited size of the simulated samples, used to determine the selection 
and particle identihcation efficiencies, is considered as systematic uncertainty. In addition, 
the hardware trigger is a potential source of systematic uncertainty due to imperfections in 
the description of data by simulation. A data sample of D*^ —)■ D^ir^, with )■ K~tt^, 

decays is used to characterise the trigger efficiencies of the pions and kaons, separated 
according to particle charge, as a function of the transverse energy of the associated cluster 
in the hadron calorimeter 25,32 . These data-driven calibration curves are used to weight 


simulated events in order to determine the efficiency of the hadron trigger. 

The effective lifetimes of meson reconstructed in a particular decay depend on 
the CP-admixture of the hnal state because CP-even and CP-odd eigenstates may have 
different lifetimes [^. Since the selection efficiency depends on decay time, this might 
lead to a source of uncertainty in the measurement. The relative change in efficiency with 
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Table 2: Systematic uncertainties on the relative branching fraction measurement for the two 
categories. The uncertainties are quoted as fractional contributions of the relative branching 
fraction and the total is the sum in quadrature of all contributions. 


Source 


B{B'^^K^K*'^) 


r+TT ) 

B{B0^K^i 

T+TT ) 

Downstream Long Downstream Long 

Fit 

0.05 

0.03 

0.20 

0.28 

Selection efficiency 

0.08 

0.10 

0.08 

0.11 

PID efficiency 

0.01 

0.01 

0.01 

0.01 

Trigger 

0.07 

0.07 

0.02 

0.09 

Lifetime 

0.05 

0.05 

- 

- 

Total 

0.13 

0.14 

0.22 

0.31 

fs/fd 

0.06 

0.06 

- 

- 


respect to the nominal value, estimated for the extreme ranges of possible effective lifetime 
distributions, is assigned as the systematic uncertainty. 

Finally, the uncertainty from the measurement of the fragmentation fractions ratio, 
fs/fd [M-16 , is taken into account. A summary of the relative uncertainties on the ratio 


of branching fractions is given in Table The final results reported in Sec. [^take into 
account correlations between the two samples; thus the systematic uncertainty for the 
combined measurement is reduced. 


6 Summary and conclusion 


A search for —)• K^K*^ decays is performed by the LHCb experiment using pp data 

recorded at a centre-of-mass energy of 7 TeV, corresponding to an integrated luminosity 
of 1.0 fb~^. The branching ratios are determined using the —)■ Arg7r’''7r“ decay as a 
normalisation mode. The measurements are performed separately for the downstream and 
long ATg categories and then combined following Refs. 

The Bg decay is observed for the first time, with a total significance of 7.1 standard 
deviations. The relative branching fraction is 


34 35 


B{B^^ ATOvr+TT-) 


0.33 ± 0.07 (stat) ± 0.04 (syst) ± 0.02 {fs/fd)- 


For the B^ decay, an upper limit at 90% (95%) confidence level (CL) is determined. 
The likelihood function is convolved with a Gaussian function with standard deviation 
equal to the total systematic uncertainty, and the upper limit is taken to be the value 
of the relative branching fraction below which 90% (95%) of the total integral of the 
likelihood function over non-negative branching ratio values is found. The central value 
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and the upper limit on the relative branching fraction of the decay —)■ K^K*^ are 


B{B^^ K^^n+n-) 

< 


0.005 ± 0.007 (stat) ± 0.001 (syst), 
0.020 (0.021) at 90% (95%) CL. 


The absolute branching fractions, calculated using the reference value of B{B^ — 
/L°7r+7r“) = (4.96 ± 0.20) x 10“^ 36 , determined without using the correlated LHCb 


measurement. The results are expressed in terms of the sum of hnal states containing 
either or mesons 


B{Bl iL°iL*(892)°) + i3(S° ^ iL°iL*(892)°) = (16.4 ± 3.4 ± 1.9 ± 1.0 ± 0.7) x 10“®, 

B{B^ ^ K^K*{m2f) + B{B^ ^ K^K*{m2f) = (0.25 ± 0.34 ± 0.05 ± 0.01) x 10"®, 

< 0.96 (1.04) X 10"® at 90% (95%) CL, 

where the first uncertainty is statistical, the second systematic, the third due to the ratio 
of the fragmentation fractions and the fourth due to the uncertainty on the branching 
fraction of the normalisation decay. These results are in agreement with theoretical 
predictions [6]-[^ and can be used to further constrain phenomenological models. 
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